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We have performed both theoretical and experimental study of optical response of parent per- 
ovskite manganites RMnOs with a main goal to elucidate nature of clearly visible optical features. 
Starting with a simple cluster model approach we addressed the both one-center (p-d) and two-center 
(d-d) charge transfer (CT) transitions, their polarization properties, the role played by structural 
parameters, orbital mixing, and spin degree of freedom. Optical complex dielectric function of single 
crystalline samples of RMn03 (R=La, Pr, Nd, Sm, Eu) was measured by ellipsometric technique at 
room temperature in the spectral range from 1.0 to 5.0 eV for two light polarizations: E || c and 
E 1 c. The comparative analysis of the spectral behavior of ei and £2 is believed to provide a more 
reliable assignment of spectral features. We have found an overall agreement between experimental 
spectra and theoretical predictions based on the theory of one-center p-d CT transitions and inter- 
site d-d CT transitions. Our experimental data and theoretical analysis evidence a dual nature of 
the dielectric gap in nominally stoichiometric matrix of perovskite manganites RMnOs, it is formed 
by a superposition of forbidden or weak dipole allowed p-d CT transitions and inter-site d-d CT 
transitions. In fact, the parent perovskite manganites RMnOs should rather be sorted neither into 
the CT insulator nor the Mott-Hubbard insulator in the Zaanen, Sawatzky, Allen scheme. 



I. INTRODUCTION 

The nature of the low-energy optical electron-hole ex- 
citations in the insulating transition metal id oxides rep- 
resents one of the most important challenging issues for 
these strongly correlated systems. All these excitations 
are especially interesting because they could play a cen- 
tral role in multiband Hubbard models used to describe 
both the insulating state and the unconventional states 
developed under electron or hole doping. Because of the 
matrix clement effect the optical response does provide 
only an indirect information about the density of states. 
Nevertheless it remains one of the most efficient tech- 
nique to inspect the electronic structure and energy spec- 
trum. 

It is now believed that the most intensive low-energy 
electron- hole excitations in insulating 3d oxides corre- 
spond to the charge transfer (CT) transitions while differ- 
ent phonon-assisted crystal field transitions are generally 
much weaker. Namely the CT transitions are considered 
as a likely source of the optical and magnetooptical re- 
sponse of the 3d metal-based oxide compounds in a wide 
spectral range 1-10 eV, in particular, the fundamental 
absorption edge. The low-energy dipole-forbidden d-d 
orbital excitations (or crystal field transitions) are char- 
acterized by the oscillator strengths which are smaller by 
a factor 10 2 — 10 3 than that for the dipole-allowed p-d 
CT transitions and usually correspond to contributions 
to the dielectric function £2 of the order of 0.001-0.01. 

Despite CT transitions are well established concept in 
the solid state physics, their theoretical treatment re- 
mains rather naive and did hardly progress during last 
decades. Usually it is based on the one-electron approach 
with some 2p-3d or, at best, 2p— > 3dt2 g , 2p— >3de g CT 



transitions in 3d oxides. In terms of the Hubbard model, 
this is a CT transition from the nonbonding oxygen band 
to the upper Hubbard band. But such a simplified ap- 
proach to CT states and transitions in many cases ap- 
pears to be absolutely insufficient and misleading even 
for qualitatively explaining the observed optical and mag- 
netooptical properties. First, one should generalize the 
concept of CT transitions taking into account the conven- 
tional transition between the lower and upper Hubbard 
bands which corresponds to an inter-site d-d CT transi- 
tion (intersite transition across the Mott gap). 

Several important problems are hardly addressed in 
the current analysis of optical spectra, including the rel- 
ative role of different initial and final orbital states and 
respective CT channels, strong intra-atomic correlations, 
effects of strong electron and lattice relaxation for CT 
states, the transition matrix elements, or transition prob- 
abilities, probable change in crystal fields and correlation 
parameters accompanying the charge transfer. 

One of the central issues in the analysis of electron- 
hole excitations is whether low-lying states are comprised 
of free charge carriers or excitons. A conventional ap- 
proach implies that if the Coulomb interaction is effec- 
tively screened and weak, then the electrons and holes are 
only weakly bound and move essentially independently as 
free charge-carriers. However, if the Coulomb interaction 
between electron and hole is strong, excitons are believed 
to form, i.e. bound particle-hole pairs with strong corre- 
lation of their mutual motion. 

Our work was stimulated by the lack of detailed and re- 
liable studies of electron- hole excitations and of a proper 
understanding of the relative role of p-d and d-d CT 
transitions in rare-earth perovskite manganites RM11O3 
(R= rare earth or yttrium) to be parent systems for the 
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colossal magnetoresistive— and multiferroic— materials. 
Given the complex phase diagram of this class of materi- 
als, studies of the nominally stoichiometric parent com- 
pound could give an insight into physics governing the 
doped version of these manganite oxides. 

Virtually all the experimental optical data available 
for manganites are focused on the LaMn.03. Except for 
paper by Kim et al£ on the absorption spectroscopy of 
thin films of RMn0 3 (R=La, Pr, Nd, Gd, Tb), there 
has been no effort to perform a comparative study of the 
optical response for different rare-earth perovskite man- 
ganites. The optical conductivity spectrum of LaMn03 
exhibits two broad intensive bands centered around 2.0 
and 4-5 eV.— ~— However, it has remained unclear just 
what the nature of the related electron-hole excita- 
tions. Some authors^i£i£ assign these both features to 
the dipole-allowed p-d CT transitions like t\ g e g — t\ g e 2 g L 
and t\ g e l g — ^2g e <?— (— denoting a ligand hole), respec- 
tively. However, others^ assign the low-energy peak to 
the Jahn- Teller orbiton-excitation 5 E g — 5 E g transition, 
or doubly-forbidden (parity and orbital quasimomen- 
tum) d-d-like crystal-field transition between two 5 E g - 
sublevels separated by a splitting due to a low-symmetry 
crystalline field. Both interpretations being particularly 
qualitative suffer from many shortcomings and give rise 
to many questions concerning the details of the charge 
transfer states or expected extremely weak intensity for 
the d-d crystal field transitions. Pronounced tempera- 
ture rearrangement of the optical spectral weight both for 
low- and high-energy bands was uncovered by Quijada et 
aL— . These authors were seemingly the first who made 
a valid conclusion that the dominant contribution to the 
optical spectral weight of the conductivity peak at 2.0 eV 
is provided by the CT hopping between nearest-neighbor 
manganese ions, or inter-site d-d CT transitions. 

Until recently main body of the optical data for man- 
ganites was obtained by reflectivity measurements fol- 
lowed by a Kramers-Kronig (KK) transformation usually 
accompanied by a number of unavoidable uncertainties as 
regards the peak positions and intensities of weak spec- 
tral features. It often implies a parasitic contribution due 
to a deterioration of the sample surface^, that can give 
rise to some ambiguities due to problems with KK trans- 
formation. The technique of ellipsometry provides signif- 
icant advantages over conventional reflection methods in 
that it is self-normalizing and does not require reference 
measurements, and optical complex dielectric function 
e = e± — ie2 is obtained directly without a Kramers- 
Kronig transformation. The comparative analysis of the 
spectral behavior of e\ and £2 is believed to provide a 
more reliable assignement of spectral features. 

First ellipsometric measurements for the single crys- 
talline LaMnOa samples were performed by Loshkareva 
et al*£ at room temperature and for the spectral range 
1.0-5.0 eV. Later on the ellipsometric measurements were 
performed for an untwinned crystal of LaMn03 by Koval- 
eva et al. n i 12 in a wide temperature range. The authors 
have presented a detailed quantitative analysis of the 



pronounced redistribution of the spectral weight (SW) 
near the Neel temperature. They concluded that the 
low-energy optical band around 2 eV consists of three 
distinct bands all assigned to intersite d-d CT transi- 
tions, and that LaMn03 is a Mott-Hubbard rather than 
a charge transfer p-d insulator as argued earlier (see 
e.g., Refsfl3lfl4T). A similar interpretation of spectral fea- 
tures near 2 eV observed in multiferroic TbMn03 was 
reported very recently by Bastjan et al.—. Many re- 
searchers— pointed to a fine sub-peak structure 
of the low-energy 2 eV band. Early optical transmission 
spectra of the LaMnOs films— revealed a fine structure 
with two features near 1.7 and 2.4 eV, respectively, which 
were assigned to the Mn 3+ d-d crystal-field transition 
5 E g — 3 Ti g , split by the JT effect. Up to this point, there 
has been little effort to understand this mysterious sub- 
peak structure. First theoretical cluster model analysis 
of the p-d CT transitions in the Mn06 octahedral did 
predict several low-energy forbidden or relatively weak 
allowed electro-dipole transitions with an onset at 1.7 eV 
to be a "precursor" of a strong dipole-allowed p-d CT 
transition at 4.5-4.7 eV. 

Such an ambiguity leaves the question of the nature of 
the main optical transitions and low-lying electron states 
in LaMn03 far from being resolved. Furthermore, it 
seems that there is a missing qualitative aspect of the 
problem that so far escaped identification. 

The one-electron band models, including such modern 
modifications as LDA+/7, do not elucidate the matter, 
first of all since they fail to reproduce important intra- 
atomic correlation effects which form the term structure 
both for the ground and excited CT configurations . 18 i 19 
In this connection we would like to emphasize the priority 
of the cluster model of the CT transitions which showed 
itself to advantage in explaining the optical and mag- 
netooptical spectra of orthoferrites, iron garnets, and a 
number of dielectric cuprates and manganites (see, e.g., 
RefsElHO-El). The model is marked by the physical 
clearness, possibility of detailed account for the electron 
correlation and the crystal field effects as well. 

Our efforts were focused on the theoretical and exper- 
imental studies of the low energy CT bands peaked near 
2 eV with a particular interest in the fine-structure ef- 
fects. To this end we first address a theoretical analysis 
of a large variety of different CT states and CT transi- 
tions in the perovskite manganites based on the distorted 
MnOe octahedra being a basic element of their crystalline 
and electronic structure. Then we present the results of 
optical ellipsometry for a number of manganites RMn03 
(R= Pr, Nd, Sm, Eu) with a supposedly reduced spectral 
weight of the d-d CT transitions and a better manifesta- 
tion of the fine-structure effects. 

The rest of the paper is organized as follows. In Sec. 
II we shortly address the electronic structure, the energy 
spectrum of a MnOe cluster, and the one-center p-d CT 
transitions. Two-center d-d CT transitions are consid- 
ered in Sec. III. The experimental results of ellipsometric 
optical measurements for perovskite manganites RM11O3 
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(R=La, Pr, Nd, Sm, Eu) are presented in Sec. IV. In 
Sec. V we address an extended discussion of experimental 
data which includes a Lorentzian fitting of dielectric func- 
tion, and a short overview of electronic structure for par- 
ent perovskite manganites as derived from optical data. 



II. ONE-CENTER (p-d) CHARGE TRANSFER 
TRANSITIONS IN ORTHORHOMBIC 
MANGANITES 



A. Electronic structure of octahedral 
3e?-metal-oxygen MeOs centers in perovskites 

The electronic states in strongly correlated 3 c? oxides 
manifest both significant correlations and dispersional 
features. The dilemma posed by such a combination is 
the overwhelming number of configurations which must 
be considered in treating strong correlations in a truly 
bulk system. One strategy to deal with this dilemma is 
to restrict oneself to small clusters, creating model Hamil- 
tonians whose spectra may reasonably well represent the 
energy and dispersion of the important excitations of the 
full problem. Naturally, such an approach has a num- 
ber of principal shortcomings, including the boundary 
conditions, the breaking of local symmetry of boundary 
atoms, and so on. Nevertheless, this method provides 
a clear physical picture of the complex electronic struc- 
ture and the energy spectrum, as well as the possibility 
of quantitative modelling. In a certain sense the clus- 
ter calculations might provide a better description of the 
overall electronic structure of insulating 3d oxides than 
the band structure calculations, mainly due to a better 
account for correlation effects. 

Five Me 3 c? and eighteen oxygen O 2p atomic orbitals 
in octahedral MeOg complex with the point symmetry 
group Oh form both hybrid Me 3g?-0 2p bonding and an- 
tibonding e g and t 2g molecular orbitals (MO), and purely 
oxygen nonbonding a lg (a), ti g (w), t lu (a), ti u (n), t 2u {^) 
orbitals. Nonbonding t\ u (o~) and ti M (7r) orbitals with the 
same symmetry are hybridized due to the oxygen-oxygen 
O 2p7r - O 2p7r transfer. The relative energy position of 
different nonbonding oxygen orbitals is of primary im- 
portance for the spectroscopy of the oxygen-3d-metal 
charge transfer. This is firstly determined by the bare 
energy separation Ae2 P 7r<T = £2p-n ~ t2p<r between O 2pit 
and O 2pa electrons. Since the O 2pa orbital points to- 
wards the two neighboring positive 3c? ions, an electron in 
this orbital has its energy lowered by the Madelung po- 
tential as compared with the O 2pn orbitals, which are 
oriented perpendicular to the respective 3c?-0-3d axes. 
Thus, Coulomb arguments favor the positive sign of the 
7r — (7 separation e p7! — e pa which numerical value can be 
easily estimated in frames of the well-known point charge 
model, and appears to be of the order of 1.0 eV. In a first 
approximation, all the 7(77) states t\ g (ir), ti u (ir), ^(tt) 
have the same energy. However, the O 2pit-0 2pix trans- 
fer yields the energy correction to bare energies with the 
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FIG. 1: (Color online) The diagram of Me 3rf-0 2p molecu- 
lar orbitals for the MeOe octahedral center. The O 2p - Me 
3c? charge transfer transitions are shown by arrows: strong 
dipole-allowed a — a and ir — n by thick solid arrows; weak 
dipole-allowed -k — a and a — n by thin solid arrows; weak 
dipole-forbidden low-energy transitions by thin dashed ar- 
rows, respectively. 



largest value and positive sign for the ii g (7r) state. The 
energy of the ti u (ir) state drops due to a hybridization 
with the cation 4pt lu (ir) state. In other words, the £i ff (7r) 
state is believed to be the highest in energy non-bonding 
oxygen state. For illustration, in Fig.l we show the en- 
ergy spectrum of the 3c?-2p manifold in octahedral com- 
plexes like Me06 with the relative energy position of the 
levels according to the quantum chemical calculations 25 
for the FeOg - octahedral complex in a lattice environ- 
ment typical for perovskites like LaFe03 and LaMn03. 
It should be emphasized one more that the top of the 
oxygen electron band is composed of O 2pir nonbonding 
orbitals that predetermines the role of oxygen states in 
many physical properties of 3c? perovskites. 

Conventional electronic structure of octahedral MnOg 
complexes is related with the configuration of the com- 
pletely filled O 2p shells and partly filled Mn 3c? shells. 
The typical high-spin ground state configuration and 
crystalline term for Mn 3+ in octahedral crystal field or for 
the octahedral MnOg - center is t 2g e g and 5 E gi respec- 
tively. Namely this orbital doublet results in a vibronic 
coupling and Jahn- Teller (JT) effect for the MnOg - cen- 
ters, and cooperative JT ordering in LaMnOa. In the 
framework of a crystal field model the 5 E g term origi- 
nates from the (3c? 4 ) 5 D term of free Mn 3+ ion. 

Unconventional electronic configuration of octahedral 
MnOe complexes is related with a charge transfer state 
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with one hole in the O 2p shells^. The excited CT con- 
figuration 7„ 3d n+1 arises from the transition of an 
electron from the MO predominantly anionic in nature 
(the 72 P hole in the core of the anionic MO being hereby 
produced), into an empty 3d type MO fag or e g ). The 
transition between the ground configuration and the ex- 
cited one can be presented as the 72p — > 3d (t 2g , e ff ) CT 
transition. 

The CT configuration consists of two partly filled sub- 
shells, the ligand j 2p and the cation 3d(t 7 2 l ge g 2 ) shell, 
respectively. The latter configuration in the case of CT 
states in the MnOg~ center nominally corresponds to 
Mn 2+ ion. It should be emphasized that the oxygen hole 
having occupied the non-bonding j 2p orbital does in- 
teract ferromagnetically with the 3d(^2ge™ 2 ) shell. This 
rather strong ferromagnetic coupling results in Hund 
rule for the CT configurations with the high-spin ground 
states. The maximal value of the total spin for the Hund- 
like CT state in MnOg - center equals 5 = 3, that points 
to some perspectives of unconventional magnetic signa- 
tures of these states. 



B. Electric dipole matrix elements 
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FIG. 2: (Color online) Two-site dipole matrix elements 
against Me3d-0 2p separation. The arrow near 4 a.u. points 
to typical Me-0 separations. For illustration we choose both 
relatively large integrals (3d z 2 \z\2p z ) governed by the Me3d- 
2p cr-bond and the relatively small ones (3d xz \z\2p x ) and 
(3d xz \x\2p z ) governed by the Me3rf-02p 7r-bond. We make 
use of hydrogen-like radial wave functions with the Clementi- 
Raimondi effective charges^ Z^ f 2 f p ^4A5 and Z^ d =10.53 
typical for Mn 3+ ion. 



The one-electron electric dipole matrix elements for 
MeOe octahedral center can be written with the aid of 
Wigner-Eckart theorem^ as follows (see Refll4l for de- 
tails) 

{w\d 9 \w) = {-iy^-»(j ll h ; j)*(7uiidii7 S ), 

\ r ? f / 

(1) 

where / \ is the Wigner coefficient for the cubic 

point group Oh, j(T) the so-called quasimomentum num- 
ber, (7u||d||7 s ) is the one-electron dipole moment subma- 
trix element. The 3d-2p hybrid structure of the even- 
parity molecular orbital:7 g /i = N lg (id r y g ^i + A 7g 2p7 g /i) 
and a more simple form of purely oxygen odd-parity 
molecular orbital 7„/i = 2p r y u [i both with a symmetry su- 
perposition of the ligand O 2p orbitals point to a complex 
form of the submatrix element in ([T]) to be a sum of local 
and nonlocal terms composed of the one-site and two-site 
(d-p and p-p) integrals, respectively. In the framework 
of a simple "local" approximation-^, that implies the full 
neglect all many-center integrals: 



(t 2 „(7r)||d||e fl ) = 0; (t 2u (n)\\d\\t 2g ) = 
(t lu (a)\\d\\t 2g ) =0; (t lu (a)\\d\\e g ) = 

{t lu (7r)\\d\\e g )=0; (t lu (Tr)\\d\\t 2g ) = 




Here, A CT ~ t pda /A pd , X n ~ t pd7I /A pd are effective 
covalcncy parameters for e g ,t 2g electrons, respectively, 
d = cRq is an elementary dipole moment for the cation- 
anion bond length Rq. We see, that the "local" approx- 
imation results in an additional selection rule: it for- 
bids the a — > 7T, and tt —¥ a transitions, ti u (er) — > t 2g , 
and ti t2u (ir) — > e g , respectively, though these are dipole- 
allowed. In other words, in frames of this approximation 
only er-type (t lu (a) e g ) or 7r-type (h t2u (ir) -> t 2g ) CT 
transitions are allowed. Hereafter, we make use of the 
terminology of "strong" and "weak" transitions for the 
dipole-allowed CT transitions going on the a — a, tt — tt, 
and 7T — a, a — tt channels, respectively. It should be 
emphasized that the "local" approximation, if non-zero, 
is believed to provide a leading contribution to transi- 
tion matrix elements with corrections being of the first 
order in the cation-anion overlap integral. Moreover, the 
nonlocal terms are neglected in standard Hubbard-like 
approaches. In Fig. [5] we do demonstrate the results of 
numerical calculations of several two-site dipole matrix 
elements against 3d metal - oxygen separation Rjv/eO- 
It is clearly seen that given typical cation-anion separa- 
tions Rji/ e o ~ 4 a.u. we arrive at values less than 0.1 
a.u. even for the largest two-site integral, however, their 
neglect should be made carefully. Exps.Q,© point to 
likely extremely large dipole matrix elements and oscil- 
lator strengths for strong p-d CT transitions, mounting 
to dij ~ eA and / ~ 0.1, respectively. 
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C. One-center p-d CT transitions in manganites 

Now we can apply the model theory to the undoped 
stoichiometric perovskite manganites RMnC>3. For our 
analysis to be more quantitative we make two rather 
obvious model approximations. First of all, one as- 
sumes that for the MnOg - centers in RMn03 as usually 
for cation-anion octahedra in 3d oxide o 20 i 25 i 29 the non- 
bonding ti g (Tr) oxygen orbital has the highest energy and 
forms the first electron removal oxygen state. Further- 
more, to be definite we assume that the energy spectrum 
of the non-bonding oxygen states for Mn 3+ Og~ centers 
in RMn03 coincides with that calculated in Refill for 
Fe 3+ Og~ in isostructural orthoferrite LaFeC>3, in other 
words, we have (in eV): 

A(ii s (7r) - ia„(7r)) « 0.8 ; A(t lg (ir) - h u (w)) « 1.8 ; 

A(h g (ir)-t lu {a)) «3.0. 

Secondly, we choose for the Racah parameters B and C 
the numerical values typical for Mn 2+ , 0.12 and 0.41 eV, 
respectively. The crystal- field parameter Dq = 0.15 eV 
provides a reasonable explanation of the Mn 2+ spectra 
in MnO 30 . Furthermore, the photoemission data^i are 
believed to confirm the relevance of this value for crystal 
field splitting parameter in RMnC>3. This set of param- 
eters is used for the model theoretical simulation of the 
overall p-d CT band in LaMn03. 14 Firstly, we argue that 
the lowest in energy spectral feature observed in LaMnC<3 
near 1.7 eV (see, e.g. Ref@) is believed to be related 
with the onset of the series of the dipole-forbidden p-d 
CT transitions £i ff (7r) —¥ e g ,t2 g , rather than with any d- 
d crystal field transition. The energy of this transition 
was picked out to be a starting point to assign all other 
p-d CT transitions. 

Weak dipole-allowed ir — a CT transitions t2 U (tt) ™ e s 
and £i u (7t) — e g form more intensive CT bands starting at 
higher than the preceding series energies, near 2.5 and 3.5 
eV, respectively, in accordance with the magnitude of the 
*iff( 7r ) — *2u(tt) an d tig(n) — ti u (n) separations. Actually, 
the ti u (ir) — e g transition at 3.5 eV has to be more in- 
tensive because the ti u (ir) state is partly hybridized with 
the t\ u (a) one, hence this transition borrows a portion of 
intensity from the strong dipole-allowed t\ u (a) — e g CT 
transition. 

The latter a — a transition forms intensive broad p- 
d CT t\ g e g — > t2 g e^t lu (a) band starting with the main 
5 E g - 6 A lg ; 5 T lu peak at s» 4.7 eV and ranging to the 
5 E g — 4 A2 g ; 5 T2u peak at w 10.2 eV with interstitial peaks 
at 8.0 eV being the result of the superposition of two 
transitions 5 E g - 4 A lg ; 5 T lu and 5 E g - 4 E g ; 5 T U , and at 
~ 8.8 eV due to another 5 E g — 4 E g ; 5 T u transition, re- 
spectively. Thus, the overall width of the p-d CT bands 
with final ^ig e g configuration occupies a spectral range 
from 1.7 up to ~ 10 eV. 

Strong dipole-allowed ir — ir CT transitions 
t2u (tt); ti u (n) — t2 g form two manifolds of equally 



intensive t\ g e g ^2g e gii,2u( n ) CT bands shifted with 
respect each other by the t2«('7r) — ti u (n) separation 
(w 1.0 eV). In turn, each manifold consists of two 
triplets of weakly split and equally intensive p-d CT 
bands related with 5 E g - 4 T lg ; 5 T U and 5 E g - 4 T 2g ; 5 T u 
transitions, respectively. In accordance with the as- 
signment of crystal-field transitions^ in LaMn03 (see 
Fig. 1) we should expect the low-energy edge of the 
dipole-allowed ir — ir CT band starting from m 4.5 eV 
(1.7 + 2.0+ (tig (ir) —t 2u ('n') separation)). Taking account 
of strong configuration interaction we should expect the 
high-energy edge of this band related with the highest 
in energy 4 T2 g term of the 3 c? 5 configuration to be 
situated near w 9.9 eV. In between, in accordance with 
our scheme of energy levels we predict peaks at 5.2, 5.5, 
6.2 (x2), 7.2(x2), 7.9, 8.2, 8.3, and 8.9 eV. The weak 
dipole-allowed a — it transitions occupy the high-energy 
spectral range from 6.7 to 11.1 eV. 

Overall, our analysis shows the multi-band structure of 
the p-d CT optical response in LaMn03 with the weak 
low-energy edge at 1.7 eV, related with the forbidden 
^ig( 7r ) — e g transition and a series of intensive bands in 
the range 4.6 + 10.2 eV starting with a composite peak 
at ~ 4.5 -7- 4.7 eV and closing by a composite peak at 
8-t- 10 eV both resulting from the superposition of strong 
dipole-allowed 7r — n and a — a CT transitions. 

In conclusion one should emphasize several points. 
First, the p-d CT excitons can move over the lattice thus 
forming an excitonic band. Its width would strongly de- 
pend on the character of the O 2p hole state. We antici- 
pate rather narrow bands for 7r — n and ir — a excitons, 
but rather wide bands for a — a and a — 7r excitons. It 
is worth noting that besides the charge transfer the p-d 
CT transitions in MnOe octahedra are characterized by 
a remarkable effect of the transfer of the spin density and 
orbital degeneracy from Mn ion to surrounding oxygen 
ions. In addition, we would like to underline a crucial role 
of the intra-atomic correlations in forming the optical re- 
sponse in 3d oxides. Indeed, the one-electron energy level 
scheme (see, Fig. 1) points to low-lying 3dt2 g orbital to 
be a final state for any lower-energy CT transition, while 
actually the lower-energy CT transitions are related with 
3de g orbital as a final state because of a prevailing gain 
in the Hund energy. 



D. Role of the Jahn- Teller distortions for the 
Mn 3+ 06 octahedra and the light polarization effects 

Orthorhombic manganites RMn03 are typical com- 
pounds with the 3c?-electron orbital ordered ground state 
where 3rf3 :E 2_ r 2-like and 3G?3 y 2_ r .2-hke orbitals for e g - 
electrons are alternately ordered on Mn 3+ sites in the ab 
plane and are stacked parallel along the c axis below the 
orbital ordering temperature Tjt ~750-1500K 32 . Such 
an orbital ordering is stabilized by the cooperative Jahn- 
Teller (JT) effect that lifts the orbital degeneracy of the 
ground state crystal 5 E term for the octahedral t^e^ elec- 
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tronic configuration of Mn 3+ ion. The orbital ordering 
causes a layer-type (so called A-type) antiferromagnetic 
AF ordering below Tjv, in which magnetic moments on 
Mn sites are aligned ferromagnetically in the ab plane 
and are stacked antifcrromagnetically along the c axis. 

The e g -electron in the JT-distorted Mn 3+ C>6 octahe- 
dral cluster occupies the superposition state 

*s = cos^|d z 2) +sin®\d x 2_ y 2) , (3) 

where the orbital mixing angle is determined by the 
deformations of Mn 3+ 6 octahedral cluster as follows: 



tan 9 = 



V3(l X — ly) 

2l z lx ly 



(4) 



where l x ,y,z are the length of the Me-0 bonds along the 
respective local co-ordinates. The unoccupied, or e g -hole 
state of the same Mn 3+ 06 octahedral cluster is described 
by the wave function ^ e which is orthogonal to ^ g : 



e 



e 



sin — 14 2 ) +cos—\d x 2_ y 2 



(5) 



and I is a mean Mn-0 bond length. For rare-earth or- 
thorhombic manganites the Mn-O/ bond corresponds to 
the medium length while two Mn-Oy bonds correspond 
to long and short lengths (see Refl32l and Fig. 3). Inter- 
estingly that if the local z-axis is directed along Mn-O/ 
bond and re-axis along Mn-O/j bond with the longest 
length, (l x — I) ^> (l z — I) fa (l y — I) and tan0 ps \/3, 
or6si 120°. In other words, it means that ^> g « \d x 2). 
In the tetragonal limit ((l x — I) (l z — t) = (l y — I)) 
the e g level splits to two singlets (Id,^) and \d y 2_ z 2), re- 
spectively), while the ti,t 2 levels split to singlets and 
doublets (see Fig. 1). Thus all the one-center p-d CT 
transitions £1,2 — > e g from filled nonbonding oxygen or- 
bitals to empty e g orbital, both allowed and forbidden, 
are believed to reveal a singlet- (quasi) doublet structure. 
The intensity of weakly dipole allowed £2« (tt) — > &g tran- 
sitions centered near 2.5 eV, can be enhanced due to 
a t2u(ft) — ti u (ir),ti u (a) mixing resulting from the low- 
symmetry deformations of MnC>6 octahedra. It should be 
noted that the probability for the £^2 — > t 2g transitions 
does not depend on the orbital mixing angle 0. 



The Jahn- Teller splitting energy is determined as follows: 
Ajt = 2gpo: where g is a vibronic constant, 



Po = [2(l x ~ l y ) 2 + 6(l z - I) 2 



(6) 



Let address the polarization properties of the main 
dipole-allowed one-center p-d CT transition t\ u (a) — > e g . 
To this end we find the dipole transition matrix elements 
in the local coordinates of MnOg octahedron: 



(tiux(cr)\d x \d z 2) = (ti U y(a)\d y \d z 2) = --(h uz (a)\d z \d z 2) 



and 



--^=-; (tx ux (a)\d x \d x 2_ y 2) = -(t luy (a)\dy\d x 2_ y 2} = -^=- . 

(7) 



/ / \ 1 7 it \ \(jd . . O r— @ . , , . 1 1 . \(jd , ("") r— O . , , . . . . QXfjd , O 

{tiux{cr)\d x \W e ) = — ^=-(sin— + V3cos— ); (t luy (a)\dy\^ e ) = — |(sin — -V3cos — ); (£i« 2 ((r)|d z |W e ) = — -^=- sin — . 

(8) 

Hence for spectral weights we obtain the following relation: 

SW X :SW y :SW Z = (2 + 3sin0 + cos©) : (2 - 3sin0 + cos0) : 2(1 - cos0) . (9) 



r 



Neglecting the tilting of MnOg octahedra we arrive at 
the ratio of spectral weights SW C and SW„t for light 
polarized || c-axis and || a6-plane, respectively: 



SW c /SW ab 



2(1 - cos©) 
(2 + cos0) ' 



(10) 



This relation was obtained earlier by Tobe et a/. 13 though 
the authors tried to apply it to explain the 2 eV band po- 
larization properties. Near the expected value of the or- 
bital mixing angles (0 > 108°) SW c /SW ab > 1.5, while 
the account for the tilting of MnOg octahedra is believed 
to reduce the ratio. 



Generally speaking, the CT exciton creation is accom- 
panied by an excitation of lattice modes. Indeed, the 
electron transfer from O 2p to Mn 3c? state, or from an 
ionic to a covalent configuration such as a — u transi- 
tion 5 E g — 6 Ai g ; 5 Ti u is accompanied by a significant 
shortening of the equilibrium Mn-0 bond length and 
a remarkable effect of the transfer of orbital degener- 
acy from Mn ion to surrounding oxygen ions accompa- 
nied by a strong change in the electron-lattice coupling. 
Thus, the incident photon can create a self-trapped ex- 
citon which can reemit a photon, returning either to the 
ground state or to various one-phonon or multiphonon 
excited states. The exciton-phonon interaction strongly 
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FIG. 3: (Color online) Scheme of Mni-Ojj-Mn2 bonding in 
the afe-plane of perovskite manganites 

affects the line-shape of absorption leading to Franck- 
Condon multiphonon sidebands and results in a phonon 
Raman scattering. The measurement of the Raman in- 
tensity as a function of excitation light energy is a very 
informative probe of the origin of electronic transitions. 

E. Temperature dependence of the spectral weight 
for the p-d CT bands 

Addressing the temperature dependence of the spec- 
tral weight for the p-d CT transitions we will first con- 
cern the effects of a magnetic ordering. At first sight 
for small one-center excitons we have a rather conven- 
tional transition with conservation of spin SAls-state and 
a spin density fluctuation localized inside the MnOg clus- 
ter. However, the redistribution of the spin density from 
manganese atom to the oxygen ones after the p-d CT 
transition switches on a strong ferromagnetic Heisenberg 
O 2p-Mn 3d exchange that leads to a number of temper- 
ature anomalies near the Neel temperature. At first, one 
has to expect a blue shift effect for the transition energy 
with the lowering the temperature near and below Tjv. 
Indeed, at T>Tjv, the average molecular field for the 
MnOg center turns into zero while the 3D antiferromag- 
nctic ordering is accompanied by a rise of the exchange 
molecular fields and respective spin splittings. Due to 
an order of magnitude bigger value of the O 2p-Mn 3d 
exchange as compared with the Mn 3d-Mn 3 c? exchange 



this is accompanied by an increase of the transition en- 
ergy with a maximal value of the blue shift as large as 
several tenths of eV. Additionally, one has to expect a 
strong (of the same order of magnitude) broadening of 
the excitonic line with the increase of the temperature 
due to strong fluctuations of molecular fields. It is worth 
noting that at variance with the inter-site d-d CT transi- 
tion the p-d CT transition is not accompanied by strong 
two-magnon Raman processes, that could be used for its 
identification. 

The temperature dependence of the phonon-assisted 
p-d CT transitions, such as fi s (7r) — > e g ,t2 g , normally 
forbidden by parity considerations, is usually described 
by the functional form predicted for such a type of pro- 
cess: 

SW(T)=SW(0)coth^, 

where SW(0), the spectral weight at K, embodies the 
effect of parity mixing introduced by lattice vibrations 
in the states connected by the transition, and hu! p is 
a phonon energy. Depending on the energy of the ac- 
tive od<i-parity phonon mode(s) such a mechanism can 
provide up to a twofold rise of the spectral weight with 
the temperature rise from helium to room temperatures. 
Obviously, this transition borrows a portion of intensity 
from the nearest in the energy dipole-allowed CT transi- 
tions. 

It should be noted that equally with the band-structure 
effects the electron-phonon interaction governs the line- 
shape and the temperature dependence of all the CT 
transitions, both dipole-allowed and forbidden ones. 



III. TWO-CENTER (d-d) CT TRANSITIONS IN 
MANGANITES 

A. Simple two-center CT model 

For illustration of covalence and correlation effects in 
a two-center MeO„-MeO„ molecular cluster, let first ad- 
dress its simplified version, namely a two-electron two- 
center (A-B) system with a single on-site electron state. 
Let a simple Hubbard-like Hamiltonian 



H 



'Ba 



r)+U Y^t n Ac 



riA- 



n B(T n B _ a ) 



V^^Ac 



(11) 



incorporates a spin- independent single-particle transfer V = 0; t ^ 0) the conventional single-particle MO-LCAO 
(tAB = tsA = t) and the on-site, as well as inter-site 
Coulomb repulsion. For a purely covalent bonding (U — 
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AA,BB 



AA+BB 



A=U-V 



AA-BB 



AB 



AB 



AA+BB 



AA-BB 



AB 



No CT corrections 



Single particle 
CT corrections 



Two particle 
CT corrections 



FIG. 4: (Color online) Energy level scheme for simple two- 
center system with the step-by-step inclusion of single-particle 
and two-particle CT corrections. Bold arrow point to an al- 
lowed electro-dipole inter-site CT transition, dashed arrow 
point to a forbidden inter-site CT transition 



bonding-antibonding basis set 



1 



$±(r) = ^=[<Mr)±0 B (r)] 



(12) 



provides a diagonalization of the transfer Hamiltonian, 
and we have simple symmetrized 

two-particle basis set for spin singlets with energies 
— \t\, 0, +\t\, respectively, and antisymmetrized 
function for a spin triplet. All purely covalent states 
describe a bond-centered charge ordering and are char- 
acterized by the equal mean on-site electron density dis- 
tribution, definite parity and zero mean dipole moment. 
At the same time, there are strong charge and dipole 
moment fluctuations. 

Depending on the properties of the on-site functions 
the single-particle bonding $+(r) and antibonding <fr-(r) 
MO-LCAO orbitals not only possess definite parity with 
regard a center of inversion in the two-center cluster, but 
resemble simple s- and p- orbitals, respectively. This im- 
plies an occurrence of the specific for " covalent" systems 
mechanism of the allowed electric-dipole transitions with 
the dipole matrix elements determined only by the AB 
separation 

(s\d\p) = ($ + (r)|d|$_(r))=<?R Afl 

and observed in the "longitudinal" polarization E | 
R-ab- In the opposite case of the purely ionic bonding 
(U 7^V r ^0;t = 0)we apply a conventional Heitler- 
London approach forming two one-site symmetrized con- 
figurations ^ aa = {4>a4>a\ and $bb = {4>b4>b} for spin 
singlets which energy equals U, and two two-site configu- 
rations, described by symmetrized even ^>ab — {0a</>b} 
and antisymmetrized odd ^[ab] = [4>a4>b] functions for 
spin singlet and triplet, respectively, which energies both 
equal V. All initial purely ionic states describe a site- 
centered charge ordering with a strictly definite on-site 



electron density distribution. Most interestingly that the 
^>aa and states interact due to a somewhat like an 
exchange resonance reaction AA O BB with a proba- 
bility amplitude Tab = aa\H\^> bb) , or an effective 
two-particle transfer integral. For a pair of identical cen- 
ters we can compose even- and odd-parity superpositions 
-^=(^aa ±*bb) with S- and P-like behaviour, and the 
energy splitting due to the two-particle transfer (±Tab), 
respectively. 

Assuming U > V, and \t\ < A = U — V we switch on 
a weak single-particle inter-center transfer and simulate 
a weak covalent situation. In frames of spin singlets we 
obtain three two-electron orbital states (see Fig. 4): 

Even-parity ground state: 

1 t 2 
^ 9 GS = -^=sm6(^AA+^BB)+cos9^AB; (E 9 GS « V--); 

(13) 

Odd-parity excited state: 

*%s = ^(*^ - {Els = U -T ab )- (14) 

Even-parity excited state: 
1 



^ES 



V2 



cos9(^aa + *bb) - sin6 ^ab; 



ES 



U+-+T AB ), 



where 



tan 29 



2t 
A' 



(15) 



(16) 



Here we ignore the spin-triplet state and ferromagnetic 
Heisenberg exchange coupling. 

Even-parity ground state ^ G g and odd-parity excited 
state ^~% s are coupled by a dipole-allowed charge transfer 
transition at the energy Tiuj 
element 



A + %r with a dipole matrix 



<*£sl<«s> 



2eR.AB sin ( 



(17) 



We see that the oscillator strength or spectral weight for 
such a two-center d-d CT transition can be easily related 
to the kinetic contribution Jkin — t 2 /A to the exchange 
integral in the AB pah:SW oc J kin / Ai&22. 

Even- and odd-parity excited states ^ GS and are 
also coupled by a dipole-allowed transition at the energy 
fiui = + 2 1 Tab | with a dipole matrix element 



d|# 



ESI 



2eR.AB cos 9 



(18) 



The d-d CT excitons can move over the lattice thus 
forming an excitonic band. Its width would strongly 
depend on the character of the electron-hole state. In 
particular, we anticipate rather wide bands for de fl -de g 
excitons, but rather narrow bands for dt2g-dt2 g excitons. 
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B. Classification of two-center (d-d) CT transitions 
in manganites 

Simple model approach can be easily generalized to 
describe inter-center charge transfer transitions between 
two MeO„ clusters centered at neighboring sites A and B, 
which define two-center CT excitons in 3d oxides. These 
two-center excitons may be addressed as quanta of the 
disproportionation reaction 

MeO v n + MeO^ -> MeO^" 1 + MeO v n +1 , (19) 

with the creation of electron MeO^ -1 and hole MeO^ +1 
centers. Depending on the initial and final single particle 
states all the intersite d-d CT transitions may be classi- 
fied to the e g — e g , e g — tig, t 2g — e 9 , and tig — t 2g ones. 
For the 3d oxides with id cations obeying the Hund rule 
these can be divided to so-called high-spin (HS) transi- 
tions S\S 2 S — > Si ± g 52 T and low-spin (LS) transi- 
tions S1S2S —> Si — — hS, respectively. 

The inter-center d-d CT transition in manganites 

MnOg" + MnOg~ ->■ MnOg" + MnOg°~ (20) 

implies the creation of the hole MnOg - and electron 
MnOg°~ centers with electron configurations formally re- 
lated to Mn 4+ and Mn 2+ , respectively. The HS one- 
particle d-d CT transition driven by e g — e g transfer can 
be written as follows: 

MnO^ (t 3 2g ; 4 A 2g e\- 5 E g ) + MnOl (t 3 g ; 4 A 2g e\- 5 E g ) -> 

MnO|-(4; 4 ^ 2g ) +MnOl°-(t%; 4 A 2g e 2 g ; 3 A 2g : 6 A lg ) . 

(-'.) 

The transition energy A for such an anti-Jahn- Teller 
transition is defined as follows: 

A = e, ]T + U d - V dd (22) 

that points to its dependence on the structural parame- 
ters (see ©)• 

The energies of different d-d CT transitions d ni d n2 — > 
d ni ~ 1 d ri2+1 can be easily derived using the proper 
Tanabe-Sugano diagrams for the and d n2+1 con- 

figurations. For instance, the low-energy HS d(e 9 )-d(e 9 ) 
CT transition 5 E g 5 E g ->■ 4 A 2g 6 A lg fl2TJ) sets up a rather 
wide band of the LS 5 E g 5 E g -> 4 A 2g 4 T 2 d(e 9 )-d(e 9 ) 
CT transitions ( 4 r 2 = 4 A lg , 4 E g ( 4 G), 4 E g ( 4 D), A 2g ( 4 F)) 
with the energy separations (given a cubic symmetry)—: 

AE( 4 A lg , 4 E g ( 4 G)) = 10 B + 5 C ; 
AE( 4 E g ( 4 D)) = 17 B + 5C; 

AE( 4 A 2g ( 4 F)) — 22B + 7 C , 

which do not depend on the crystal field splitting pa- 
rameter Dq and do scarcely vary from one compound to 



another that makes the transitions to be important opti- 
cal signatures of the d(e 9 )-d(e 9 ) charge transfer. Making 
use for the Racah parameters the numerical values typical 
for free Mn 2+ ion: B=0.12; C=0.41 eV, we arrive at three 
LS d(e 9 )-d(e 9 ) CT bands separated from the low-energy 
HS d(e 9 )-d(e 9 ) CT band by 3.2, 4.1, and 5.5 eV, while 
from experimental optical spectra for MnO— we obtain 
slightly lower energies: 2.8, 3.8, and 5.1 eV, respectively. 

It should be noted that the intra-atomic electron- 
electron repulsion does mix the f 2 ll g ; 2Sl+1 Fie™ 2 ; 252+1 r 2 : 
2S+ip s^^gg w ith different configurations and/or differ- 
ent intermediate momenta Si. 2 ,Ti^ 2 but the same crys- 
tal terms 2S '+ :L r. This point is of a great importance for 
4 E g ( 4 G) and 4 E g ( 4 D) terms of Mn 2+ which are a result 
of the interaction of two bare terms: t 2g ; 4 A 2g e 2 ; 1 E : 4 E g 
and t\ g 2 Ee 2 ; 3 A 2g : 4 E g , respectively. It is worth noting 
that only the former t\ g ] 4 A 2g e 2 ; X E : 4 E g term is active 
in the d(e 9 )-d(e 9 ) charge transfer, which does not affect 
the t 2g subshcll. 

The wave functions for 4 E g ( 4 G) and 4 E g ( 4 D) terms 
can be written as follows: 

V( 4 E g ( 4 D)) = 

cosa^^A^'Eg : 4 E g ) +sma^( 2 E g ; 3 A 2g : 4 E g ) 

*( 4 E g ( 4 G)) = 

smaV^A^Eg : 4 E g ) ~ cosa^{ 2 E g - 3 A 2g : 4 E g ), 

(23) 

where tan 2a = 4v3, and a ~ 41°. In other words, the 
both terms are almost equally involved in the d(e 9 )-d(e 9 ) 
charge transfer under consideration. 

In the pa 3 eV gap in between the low-energy HS 
and LS d(e 9 )-d(e 9 ) CT transitions one may observe a 
relatively weak low-energy HS d(i2 9 )-d(e 9 ) CT transi- 
tion 5 E g 5 E g -> 4 T 2g 6 A lg at the energy A + lODq and 
low-energy LS d(e 9 )-d(i 29 ) CT transitions 5 E g 5 E g -t 
4 A 2g 4 Ti g , 4 T 2g anticipated at the energies A + 2.1eV and 
A + 2.4 eV, respectively—. Spectral weights for the both 
transitions are determined by t 2g — e g( e g ~ ^2 9 ) trans- 
fer integrals which are believed to be relatively small as 
compared with large e g — e g transfer integrals. 



C. Two-center CT d-d transitions:Effect of orbital 
states and Mei-0-Me2 bond geometry 

As it is seen from (jTHJ) , (fTT|) the d-d CT transition prob- 
ability amplitude is determined by the respective single- 
particle transfer integral. These integrals depend both 
on initial and final states and the Mei-0-Me2 bond ge- 
ometry. For two octahedral MeOg clusters sharing the 
common oxygen ion we arrive at following expressions 
for transfer integrals^ 
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ii 2 (e s 0;e g 0) « t ss +t aa cos9; t 12 (e g 0; e g 2) = t 12 (e g 2; e g 0) = t 12 (e g 2; e g 2) = 0; 



ti 2 (e g 0;t 2g /i) « t a7r D^ (w); t 12 (t 2g /i; e g 0) « t^Dj^u;); t 12 {t 2g nx;t 2g \± 2 ) w <W-D^ 2 (w) 



(24) 



r 



where the 3d-orbitals e g /j,(e g = d z 2,e g 2 = d x 2_ y 2) y t 2g )i 
(t 2g ± 1 = T^|(rfxz ± id yz ),t 2g 2 = d xy ) are speci- 
fied in the local co-ordinates for MeiOe and Me2C>6 
clusters with zi and 2:2 axes directed to the common 
oxygen ion; £ S s,crcr,<r7r,7ra-,7r7r are transfer parameters for 
ss, era, air, wa, irir bonds, respectively; 6 is a Mei-0-Me2 
bonding angle; D^}^ 2 (w) is the Wigner rotation matrhs^l 



with ui being the Euler angles, specifying the transfor- 
mation from the local co-ordinates for MeiOe to that for 
Me206 cluster. 

Thus we see that the probability amplitude for the 
HS # g (l) -> * e (2) CT e g -> e g transition along c-axis 
in perovskite manganite RMn03 would depend on the 
crystallographic parameters as follows: 



A c (* fl (l) -> * e (2)) = A c (* g (2) -> * e (l)) oc <*,(l)|d,|* e (2)) oc *JML±*lWl = lsine ( ^ + ^ C ° Sgc) (25) 



r 



while for the similar transition in the pair of nearest 
neighboring Mn 3+ 06 octahedral clusters located in ab- 



plane we arrive at 



A ab ( % (l) - * e (2)) oc <*,(1)I«U*.(2)> oc = _ sin6 ) ^±^M ; (26) 



A ab (V g (2) -> * e (l)) ex (* 9 (2)|d a!iI ,|* e (l)) ex 



i(* fl (2) * e (l)) 



i ( ^ +sine) (!fl±^£M. (27) 
2 V 2 A 



Hence the spectral weight for the inter-site d-d CT tran- 
sitions may be written as follows: 



Qu r 1 ■ 2 n f tss + tag COS 9 C 

o vVc oc — sin 6 

4 V A 



SW ab <x i(J+sin 2 e) 



* ss + t cos 0. 



A 



(28) 



It is worth noting that the orbital mixing angle depen- 
dence (|2"8"|) strictly coincides with that of RefflU Kim el 
al£ presented an expression for SW ab which accounts 
only for 1 — > 2 transfer (compare our Exp. (|26[ ) and 
Exp. (4) from RefQ) that immediately questions their 
main conclusions. 

The spectral weight for the two-center de g -de g CT 
transition can be easily related with the kinetic contribu- 
tion to the e g — e g exchange integral in the Mn 3+ -Mn 3+ 



pair: 



SWr. 



SW ab ot 



■ sin 2 6; 



jab ( 



A 



^(^+sin 2 9). (29) 



The both expressions do predict a sizeable suppression of 
the SW a fc and SW C in RMn0 3 for the two-center de g -de g 
CT transition with decreasing the ionic radius of R-ion, 
if to take into account proper variation of the structural 
parameters 6, 0, and Ajt oc p— . 



D. Two-center CT d-d transitions:Spin dependence 
and temperature behaviour 

After ruling out the roles of the Mn-O-Mn bond angle 
and the octahedral-site distortion in determining the d-d 
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CT transition probability, we have to consider whether dipole moment operator its matrix on the pair wave func- 
the spin degree of freedom affect the d-d CT transition. tions depends on the spin quantum numbers 
Interestingly that despite the spinless character of 



J 



(* 1 $2S 1 S 2 SMs\d q \& 1 & l S , 1 S! i S , Ms>) = (-1) S <WW/ S , y/[Si,S' 2 \ j § § ? } (*i*al4l*i*a) , (30) 



r 



for the 1 — > 2 (Si — > S 2 ) electron transfer. Here the symbol, [5] =25+1. For the partial spectral weight for 

S1S2S — > S[S 2 S transition we have 



conventional notations 



are used for spin 6j- 



SW(SiS 2 S -> S[S' 2 S) oc [S^J^K^^^^SMsldgl^^S'^SMs)^ = [S U S 

M s 



where 



Si S*2 S 
S 2 + 5 Si - \ 5 



(-1) 



S1+S2+S 



S{S + 1) - (Si ~ S 2 ){Si - S 2 - 1) 
2Si(2Si + l)(2S 2 + l)(2S 2 + 2) 



M sl | f } K*iNrf>'i*2>l 2 

(31) 
(32) 



1/2 



5*i 5*2 5* 
a 1 o 11 
02 2 1 22 



S1+S2+S 



(Si + S 2 )(Si + S 2 + 1) - S(S + 1) 
2Si(2Si + l)2S 2 (2S 2 + l) 



1/2 



Making use of the normalization relation for 6j- 
symbols 27 



(33) 



Ere c'i j Si S 2 S 
^-^j j s 2 ^ i 

S' 



= 1 



(34) 



we arrive at a spin sum rule: 

[S]- 1 J2 K*i*aSi&5M s |d ? |$i^5i^5Ms)| a = |($ 1 *2|4|*i«2)| 2 » 

For the partial spectral weight for S1S2 — ^ S[S 2 transition we have 



(35) 



r 



(36) 



where is the temperature dependent statistical weight S(S^) would be determined by a statistical average 
of SiS 2 S spin multiplet. Taking into account the above (S(S + 1)) = (S 2 ) which, in its turn, relates to the spin- 
expressions for 6j-symbols we see that the temperature 
dependence of the partial spectral weight SW(SiS 2 — > 
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spin correlation function " ^ 2 j / as f° uows 

((Si ■ S 2 )) = \ l(S(S + 1)> - Si (Si + 1) - S 2 (S 2 + 1)] 

Thus we should conclude that the partial spectral weight 
for S1S2S — > S[S 2 S transition in an isolated spin pair 



is governed by a spin-dependent prefactor containing the 
spin-spin correlation function: ■ S 2 ^. For the HS 

transition Si Si — > Si — 5 Si + \ and the LS transition 
S\Si — > Si — hSi — 5 in the pair of identical 3d ions we 
arrive at 



QW ( Q c v c 1 c , ^ (S(S + 1)) 

SW(SiSi ^Si- -Si + -) oc 2Si(2Si + 1) 



(S 1 -S 2 ))+5 1 (5 1 + 1) 



Si(2Si + l) 



(38) 



<?™<7 S [2^(2^ + 1) -(5(5 + 1))] 
SPF(SiSi -> Si - -Si - -) oc 25l(2Sl + 1) 



5 1 2 -((S 1 -S 2 )y 



5i(25i + l) 



(39) 



1.0 
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FIG. 5: Schematic temperature dependence of the spectral 
weight for the HS d-d CT transition in ferromagnetically and 
antiferromagnetically ordering pairs, respectively 

respectively. These expressions allow to obtain eas- 
ily both the low-temperature (T <C Tjy) and high- 
temperature (T 3> Jjv, ■ Sj^j — > 0) limits for the spin 

prefactor. Schematically the temperature dependence of 
the spectral weight for the HS d-d CT transition in ferro- 
magnetically and antiferromagnetically ordering Mn 3+ - 
Mn 3+ pairs (Si = 2) is shown in FigJSJ It should be 
noted that similar results were obtained earlier by Oles 
et ai— though from another starting point. 

In accordance with a spin sum rule (|35l) the sum of spin 
prefactors in the right hand sides of the expressions (|3"5)) 
and (|39[) turns into unity due to an exact compensation 
of temperature dependent terms with the spin-spin cor- 
relation function. In other words, the higher-energy LS 
bands exhibit a strictly opposite SW evolution with the 
temperature as compared with HS band, that is vary- 
ing the temperature we arrive at the SW transfer be- 
tween the HS- and LS-subbands. In the limit T — > the 
HS band in the antiferromagnetically coupled spin pair 



or LS band in the ferromagnetically coupled spin pair 
are freezed out. However, this prohibition is lifted for 
the HS band in the antiferromagnetically coupled spin 
pair, if we take account of effective staggered magnetic 
fields which affect the spin pair in the antiferromagnetic 
surroundings giving rise to a spin singlet-triplet mixing. 
Furthermore, the low-temperature SW for the HS band 
in the antiferromagnet is believed to be a quantitative 
measure of the Neel state weight in the ground antifer- 
romagnetic state. The ratio of spectral weights SW C and 
SW a b does weakly depend on the structural parameters 
such as and 9, in the high-temperature limit (T ^> T/v) 
SW C / SW a b ~ 0.5, while in the low-temperature limit 
(T Tjv) for antiferromagnetic A- type order this ratio 
turns into zero. It is important to note that the tem- 
perature variation between these two limits is controlled 
by the nearest-neighbor spin-spin correlation function, 
which varies smoothly with T, and not by the square of 
the sublattice magnetization, which of course vanishes 
above Tjy. One generally has considerable short-range 
order above any magnetic-ordering temperature, and for 
LaMnC>3 with TV = 140 K the ferromagnetic correlations 
should persist to quite high temperatures and the full 
difference in the spectral weight of the d-d CT transi- 
tion between the antiferromagnetic ground state and the 
paramagnetic state is not saturated by 300 K. In contrast 
with LaMnOa the room temperature seems to correspond 
quite well to a high-temperature limit (T 3> Tzv) for such 
manganites as SmMn.03 and EuMnC>3 with Tjv ~ 60 and 
50 K, respectively^ 2 .. In other words, the optical measure- 
ments could be used to inspect the short-range spin cor- 
relations in manganites. It is worth noting that the spin- 
spin correlation function ^ ^Sj ■ Sj^j ^ for the nearest- 
neighbor spins localized on the i-th and j-th Mn sites 
governs the effect of phonon renormalization— . This ef- 
fect was examine d 36 ! 37 for all orthorhombic RMn03 com- 
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pounds. 



E. Some features of CT transitions in strongly 
coupled corner-shared MnOg clusters 



is a normalization factor. If we neglect the overlap inte- 
grals, the probability amplitudes for configurations cou- 
pled by a single p-d transfer are defined by the ratio of 
the generalized p-d transfer integrals and the p-d transfer 
energy as follows: 



Above we have addressed a somewhat idealized sce- 
nario of p-d and d-d CT transitions which implies well 
isolated or weakly coupled MnOg clusters. Actually in 
perovskite manganites we deal with strongly coupled 
corner-shared MnOg clusters sharing a common oxygen 
ion. Strictly speaking, it means that we cannot make use 
of an oversimplified classification of the p-d and d-d CT 
transitions. Indeed, the O 2p electrons localized on the 
intermediate common oxygen ion cannot be attributed 
to one or another cluster. Formally this implies a strong 
overlap of the bare wave functions localized on the neigh- 
boring clusters. Generally speaking, this trouble can be 
overcomed with making use of the orthogonalization pro- 
cedure, e.g. by the Lowdin technique. This immediately 
results in a two-center character of the wave functions 
with the onset of two-center spin and orbital correlations. 

Many features of the corner-shared Me06 clusters re- 
lated with the p-d overlap and covalency effects can be 
illustrated by a more simple generic two-hole-three-site 
M1-O-M2 cluster model to be a natural generalization of 
the above addressed two-site A-B model. To this end we 
make use of a technique suggested in Refs i22r— to de- 
rive the expressions for the Dzyaloshinsky-Moriya cou- 
pling and spin-dependent electric polarization in copper 
oxides. 

We start with the construction of spin-singlet and spin- 
triplet wave functions for our three-center two-hole sys- 
tem taking account of the p-d hopping, on-site hole-hole 
repulsion, and crystal field effects for the bare ground 
state 101 and bare excited configurations {n} (Oil, 110, 
020, 200, 002) with different hole occupation of Mi, O, 
and M2 sites, respectively. The p-d hopping for M-0 
bond implies a conventional Hamiltonian 



M'r's _ ({nyT'S\H\{n}TS) 



H, 



pd = tpadfiPadp + h.C. . 
a p 



(40) 



where p\ creates a hole in the a state on the oxygen site, 
while dp annihilates a hole in the /3 state on the copper 
site. 

Perturbed wave functions can be written as follows 



{n}'T'S 



}FS ®{n}T'8Ml) 

(41) 



M'^{n},r' 



where &{ n }rsM are bare wave functions for {n} configu- 
ration and the summation runs both on different config- 
urations and different orbital T states; 



-1/2 



r7 rs 

V{n} 



E 

{n}'^{n},r> 



{n}rs 



(42) 



1n}rs 



In other words, 



A({n}T'S - {n}TS) 



(43) 



„200| 



--110I' •-- r ~ w fiil e the probability amplitudes, or 

hybridization parameters for configurations coupled by a 
two-step p-d transfer, such as c 2 ™ , c?jj5 , c™? are on the 
order of |-^-| 2 . 

Formally, the 101^011(110) and 101^002(200) CT 
transitions may be attributed to p-d and d-d charge 
transfer, respectively. However, the configuration inter- 
action due to the charge transfer and overlap do result in 
that both types of the CT transitions are accompanied 
by a charge redistribution all over the three centers. Fur- 
thermore, a novel type of p 2 -dd CT transitions 101— >020 
does emerge due to a mixing of all the 101, 110, 011, 020 
configurations in the initial and final states. Obviously, 
the main contributions to dipole transition matrix ele- 
ments for p-d, d-d, and p 2 -dd CT transitions are propor- 
tional to £i ||^| 2 , and ||^| 2 , respectively Accord- 
ingly, the spectral weights for the dipole allowed p-d, d-d, 
and p 2 -dd CT transitions are 

SW p . d cx (J^j ; SW d . d cx SW P 2_ M cx 

respectively. 

Interestingly, all the dipole-allowed CT transitions 
obey the AS* = spin-selection rule, however, the dipole 
transition matrix elements depend on the total spin S due 
to a spin (S) dependence of the energy denominators in 
the probability amplitudes and so in the normalization 
factors, particularly due to a contribution of 200, 020, 
and 002 two-hole on-site configurations. Indeed, for such 
configurations the spin-singlet and spin-triplet terms are 
usually separated by a large energy gap. In general, the 
spin-dependent part of the electro-dipole moment oper- 
ator acting in the basis of the three-site wave functions 
can be written as follows 



'101 1 



.1101 

'101 1 



4)20 I 



„020| 
L 110 1 



n(si . s 2 ) 



(44) 



Such a spin operator for a transition dipole moment was 
firstly introduced by Sugano et al— to explain so-called 
magnon side-bands in the optical absorption spectra of 

magnetic crystals. Simple estimates yield |II| cx f^f^) 

and |II| cx (^-^-^j for p-d and d-d CT transtions, re- 
spectively. In other words, the spin-dependent part of 
the spectral weight for both transitions appears to be of 

the same order: SW sp i n cx (2^) ■ It should be noted 
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that the spin- dependent part of the spectral weight in 
the three-site two-hole cluster can be written as follows: 

SWspin = (SWtriplet ~ SWainglet) (§1 ■ §2) , (45) 

if to take into account that 

(§i • § 2 ) 2 = ^ - ^(si • s 2 ). 

From this reasoning we may conclude that, in contrast to 
a clue point of papers by Kovaleva et al. 11 i 12 ; the both 
dipole-allowed d-d and p-d CT transtions are equally 
sensitive to the temperature-dependent spin correlations, 
only with a different relative change of the optical re- 
sponse through the onset of the magnetic order, relatively 
large for the former and relatively small for the latter. 

IV. ELLIPSOMETRIC STUDY OF OPTICAL 
RESPONSE OF RMnO s 

Theoretical analysis (see preceding section and 
Refs fTTIfT^ ) clearly shows that the spectral weight of the 
HS de g -de g CT transitions is to be sizeably suppressed 
with decreasing the ionic radius of R-ion and increasing 
structural distortions that provides an opportunity to a 
better manifestation and inspection of the relatively weak 
p-d CT transitions expectedly forming a fine structure of 
2 eV band in manganites. Indeed, a marked suppression 
of the spectral weight for the 2 eV band was observed by 
Kim et al— in the ab-plane absorption spectra of RMn03 
(R=La, Pr, Nd, Gd, and Tb) thin films. However, the 
authors focused on the suppression effect and did not 
concern the fine structure of the band. We have, on the 
contrary, studied the optical response of RMn03 (R=La, 
Pr, Nd, Sm, and Eu) with a main interest in the fine 
structure of 2 eV band and a more comprehensive as- 
signment of different p-d and d-d CT transitions. 

The perovskite manganites RM11O3 crystallize in the 
orthorhombic structure (space group Pbmn), where 
Mn 3+ 06 octahedra share corners with Mn 3+ -0-Mn 3+ 
angle varying from w 155° for LaMn03 to w 140° in 
the end of the series—. The orthorhombic RMn03 com- 
pounds with large ionic radia from La to Gd are char- 
acterized by A-type antifcrromagnetic AFM structure. 
In this structure the Mn moments are ferromagnetically 
(FM) ordered in the ab planes and antiferromagneti- 
cally coupled along the c direction. The Neel temper- 
ature Tjv decreases from 140 K for LaMn03 to 40 K 
for GdMnOs— . With increasing structural distortions, 
the magnetic structure of RMn03 compounds with small 
ionic radii changes to an incommensurate IC sinusoidal 
AFM structure (R=Tb, Dy) or to a so-called E-structure 
(R=Ho, Er). 

We have studied the optical response for single crys- 
talline samples of RMnOs (R=La, Pr, Nd, Sm, Eu) 
grown by the floating-zone technique described in Ref|42|. 
Magnetic properties of these samples are reported in 
RefslHil As a rule, the as-grown single crystals of 



RMnOs had twins, however, it was possible to choose 
some samples exhibiting a strong magnetic and optical 
anisotropy, hence these could be considered as practically 
untwinned. 

The measurements were performed with a home-built 
automatic ellipsometer on the x-ray oriented single crys- 
talline bulk samples with as-grown polished surface (~ 
3x5 mm 2 ) where the angle of incidence was 67° or 
72°. Spectroscopic ellipsometry measurements were per- 
formed on the surfaces for PrMnC>3 and NdMn03 crys- 
tals for high-symmetry orientations a and c of the optical 
axis with respect to the plane of incidence of light and 
on the ab surfaces for SmMn03 and EUM11O3 crystals 
for high-symmetry orientation a. The crystal surfaces 
(ac-plane for PrMnOs and NdMn03, and a6-plane for 
SmMnOs and EuMnOs) were polished to optical grade 
with diamond powders. 

Optical complex dielectric function e = e± — ie2 of sin- 
gle crystalline samples of RMnOs was studied at room 
temperature in the spectral range from 0.5 to 5.5 eV. 
However, our experimental technique had some spectral 
limitations and provided reliable optical data only in 
the spectral range 1.5-5 eV of strong absorption. Op- 
tical constants, refractive index n and extinction ratio 
k have been calculated with using relations reported in 
Refl45l. The technique of ellipsometry provides signifi- 
cant advantages over conventional reflection methods in 
that it is self-normalizing and does not require reference 
measurements, and optical complex dielectric function 
e = £1 — iei are obtained directly without a Kramers- 
Kronig transformation. The comparative analysis of the 
spectral behavior of E\ and £2 is believed to provide a 
more reliable assignement of spectral features. 

Fig. [5] shows the £\ and £2 room temperature spectra 
of RMn0 3 (R= Nd, Pr, Sm, Eu) for two main light 
polarizations, E _L c-axis and E || c-axis, respectively, 
together with earlier data for LaMnO^^. The opti- 
cal spectra of RMnOs are interesting in several aspects. 
First of all, as in many earlier works on LaMn03, we 
clearly see two intensive and rather broad optical fea- 
tures peaked at around 2 eV and 4.5 eV, respectively. 
Both the intensity and shape of the high-energy band 
change weakly inbetween different manganites, whereas 
the spectral weight of the broad band peaked at 2.0 eV in 
LaMn03 seems by factor 1.5 bigger than in other man- 
ganites. This result agrees with the findings by Kim et 
al£ who investigated the afr-plane absorption spectra of 
RMn0 3 (R=La,Pr,Nd,Gd, and Tb) thin films and ob- 
served a drastic suppression of the 2 eV peak with the 
decrease of the R-ion radii. 

The broad 2 eV band exhibits strong anisotropy, it is 
hardly visible at E | c-axis polarization. As expected, 
our measurements for R=Nd, Pr, Sm, Eu have uncov- 
ered many subtle features missed in earlier optical ex- 
periments on parent manganites. First, the compara- 
tive analysis of the £\ and £2 spectra (see Fig. [5]) allows 
us to unambiguously conclude that the 2 eV feature is 
composed of a single intensive and rather broad band 
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FIG. 6: (Color online) Real and imaginary part of the dielectric function e in RMnOs. The left hand panel presents the spectra 
for E _L c polarization, the right hand panel does for E || c polarization. The Lorentzian fitting data are shown by solid lines. 
A proper Lorentzian decomposition with assignment of all the peaks is shown for £2. 



peaked at 2.0 eV for LaMn03 and 2.2 eV for other man- 
ganitcs, and several relatively weak and narrow bands 
peaked nearly equally for all the manganites at 1.3, 1.9, 
2.3, 2.7 eV (e a b), respectively. Interestingly, that not 
only the energies but also the intensities of these bands 
hardly depend on the manganite type if any. This ob- 
servation immediately points to different nature of both 
features at variance with Kovaleva et alM who addressed 
the 2 eV feature in LaMn03 to consist of three distinct 
bands, peaked at 1.95, 2.35, 2.70 (s 2a b, T=300K) each 
attributed to inter-site d-d CT transitions. Weak sub- 
peaks show up a clear anisotropy. 

As regards the high energy spectral feature peaked at 
around 4.5 eV one should note its hardly visible com- 
posite structure with as minimum three rather intensive 
bands which integral spectral weight for E || c-axis po- 
larization seems to be slightly larger than for E _L c-axis 
polarization. 

Our ellipsometric data for RM11O3 samples (R=Pr, 
Nd, Sm, Eu) agree with that of Kovaleva et alJ^- for 
LaMnC>3 samples all over the spectral range with a vis- 
ible discrepancy only below 1.5 eV where ellipsometric 
technique fails to reliably reproduce the optical response 
because of a crucial role of the sample morphology, sur- 
face roughness, and depolarization effects. Namely these 
effects rather than intrinsic electronic structure are be- 
lieved to be responsible for a sizeable scatter of £206 data 
for different RMnC-3 samples (R=Pr, Nd, Sm, Eu) below 
1.5 eV. The spectral limitations affect smoothly the re- 
sults in the low absorption range, but not the positions 
and the strength of intensive absorption bands. In other 
words, these limitations are believed not to affect the 
main result of our paper that is the assignment of p-d 



and d-d CT transitions in the spectral range 1.5-5 eV. 



V. DISCUSSION 

A. Optical response of RM11O3 as an interplay of 
p-d and d-d CT transitions 

In Fig. [5] we have presented the results of a semi- 
quantitative analysis of the spectra in the energy range 
covered by our experiment based on a dispersion analysis 
of the complex dielectric function e, which was fitted by 
a set of 8 Lorentzian oscillators. To reproduce the tails of 
the high-energy bands we have introduced an additional 
optical band peaked at around 8.7 eV (e a b) or 8.2 eV (e c ). 
It is worth noting, that the Lorentzian fitting should be 
made with a great care because of several points. First 
of all, the Lorentzian function is a clear oversimplifica- 
tion for complex line-shapes resulting from an interplay 
of the electron-lattice interaction and excitonic band ef- 
fects. Second, it is worth noting the increased uncer- 
tainties both of the measurement and, accordingly, of 
the analysis close to the low- and high-energy cutoffs of 
our experiment. It should be noted that the results of 
the Lorentzian fitting for the E\ near 2 eV optical feature 
rather than £2 are quite robust to the assumptions related 
with the transitions beyond the energy range. All the os- 
cillators we made use can be sorted into three groups 
assigned to: i) an intensive and broad band peaked at 
2.2 eV, ii) weak and narrow bands peaked at 1.3, 1.9, 
2.3, 2.65 eV, forming a fine structure of 2.2 eV band, and 
iii) intensive bands peaked at 3.0, 3.9, 4.65, 5.3 eV, form- 
ing a broad spectral feature peaked at around 4-5 eV. A 
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proper Lorentzian decomposition is shown in Fig.[6j 

Both the spectral weight, polarization properties, and 
temperature behavio r 9 ' 11 ' 12 of the intensive and broad 2 
eV band in all the manganites investigated point to its 
inter-site d-d CT character, or strictly speaking, to the 
low energy HS e g — e g CT transition. Such a conclu- 
sion is strongly supported by proper quantitative esti- 
mates-iLl 2 ^. 

One of the most impressive manifestations of the d- 
d CT origin of the 2 eV band in LaMn03 is provided 
by the resonance Raman spectroscopy—. Indeed, the 
out of phase Pj2 g (Pbnm system) breathing mode— at 
611 cm -1 , which strongly modulates the intersite d-d 
charge transfer, shows a single pronounced resonance at 
around 2 eV. It is worth noting that the authors— found 
a good agreement between the experimental results and 
a theory by Allen and Perebeinos— based on the Franck- 
Condon mechanism activating multiphonon Raman scat- 
tering. However, the assignment of the 2 eV absorption 
peak to the intra-molecular transition due to the Frank- 
Condon process between the Jahn- Teller split E g levels 
seems to be questionable because it cannot explain the 
basic features of 2 eV band such as a rather strong spec- 
tral weight with a peculiar temperature dependence (see 
detailed discussion in Ref.|H). 

As it is seen from Exp. (|29|) the spectral weight of the 
2 eV d-d CT band depends on the orbital mixing angle 
0, Mn-O-Mn bond angle 8, and the d-d transfer energy 
A, which in turn depends on the JT distortions of MnC>6 
octahedra (see Exp. (|2"2"|l ). Our experimental data unam- 
biguously show that the minimal inter-site CT energy 
A ps 2.2 eV is virtually constant for R= Nd, Pr, Sm, 
and Eu. Interestingly that the red shift of the band in 
LaMn03 (Az, a = 2.0 eV) is due to a sizeably lower mag- 
nitude of the JT energy ejt oc Tjt (see RefJH) as com- 
pared with other manganites studied (Ap r ,Nd,Sm,Eu = 
2.2 eV), furthermore, this can provide a sizeable 20% 
contribution to the spectral weight enhancement, while 
the effect of O and 8 variation is restricted by ~10%. 
This result agrees with theoretical predictions based on 
the structural dependence of the SW and also supports 
the inter-site d-d CT assignment of the 2 eV band. 

Lorentzian decomposition of 2 eV feature in RMnO.3 
with assignement of all the peaks is shown in FigO The 
intensity and energy position of four narrow bands form- 
ing the fine structure of 2 eV feature agree surprizingly 
well with the theoretical predictions made earlier by one 
of the authors (see Refil4) for the forbidden ii fl (7r) —> e g 
and weakly-allowed t2 U (7r) — > & g one-center p-d CT tran- 
sitions splitted by a near tetragonal crystal field. It is 
worth noting that a close similarity of tetragonal distor- 
tions for MnOg octahedra in all the manganites under 
consideration— does explain a nearly equal position of 
relatively weak and narrow bands peaked for all the man- 
ganites at 1.3, 1.9, 2.3, 2.7eV (e a b). Weak non-tetragonal 
distortions of MnOg octahedra are seemingly visible for 
tigin) — ► e g transition. Indeed, the a6-response is well 
fitted by two Lorentzians peaked at 1.3 and 1.9 eV, while 



c-response needs two Lorentzians peaked at 1.3 and 1.7 
eV, respectively. It is worth noting that the energy split- 
ting of nonbonding oxygen (7r), ^(tt) orbitals («0.5 
and ps 0.4 eV, respectively) appears to be slightly smaller 
than that of e g orbital (Ajt ps 0.7 eV). The relative 
sharpness of all the weak p-d CT bands agrees with a lo- 
calized character of nonbonding oxygen O 2pir orbitals. 
The most intensive and broad high-energy band mod- 
eled by four Lorentzians peaked at 3.0, 3.9, 4.65, and 
5.3 eV can be unambiguously assigned to the one-center 
dipole-allowed p-d CT transitions: ti u (ir) — e g (jx — a), 
*2u(tt) — t2g(^ — tt), ti u (a) — e g (cr — a) transitions, re- 
spectively, with a hardly resolved low-symmetry splitting 
for the most intensive and broad bands. The composite 
high-energy p-d CT band is markedly revealed by the 
resonance Raman spectroscopy in LaMn03— . One ob- 
serve pronounced resonances at around 4-5 eV both for 
the A g Jahn- Teller stretching and A g out-of-phase rota- 
tion modes of the oxygen octahedra at 496 cm -1 and 
284 cm -1 , respectively, and also for the in phase sym- 
metric Bi g stretching mode at 655 cm -1 . Interestingly 
that the same Raman modes resonate also at around 2 
eV that supports our conclusion about the composite d- 
d/p-d structure of the low-energy CT band. It is worth 
noting that in LaMnOs we deal with a near degeneracy 
of the d-d CT transition and one of the p-d CT transi- 
tions at 1.9-2.0 eV. It means a strong coupling/mixing of 
two types of transitions. 

Concluding the subsection we would like to comment 
an assignment of p-d and d-d CT transitions in LaMn03 
reported in recent papers by Kovaleva et a l 1 ' ■ As it 
is clearly seen in Fig.[6jsee also Figs. 8 and 9 in Ref.[l2l) 
the optical response of LaMnOs in contrast to RM11O3 
(R=Pr, Nd, Sm, Eu) does not distinctly reveal fine struc- 
ture effects. Neverthless, the authors—^ made use of a 
Lorentzian decomposition of the complex dielectric re- 
sponse £{, and e c , measured by ellipsometric technique in 
the spectral range 1.2-6.0 eV at the temperatures var- 
ied between 20 and 300 K. They did consider only one 
high-energy strongly dipole-allowed p-d CT transition at 
4.7 eV, while other optical bands observed at around 2.0, 
2.4, 2.7, 3.8, 4.4, 4.6, 5.7 eV were assigned to different 
HS and LS inter-site d-d CT transitions. The assignment 
was based on the temperature dependence of the par- 
tial and total spectral weights, and a parameterization of 
the energy spectrum of final d 3 -d 5 configuration. How- 
ever, the argumentation is believed to be questionable 
because of several unacceptable assumptions and incon- 
sistencies. First, the authors neglect a number of one- 
center p-d CT transitions to be precursors of the most 
intensive high-energy dipole-allowed p-d CT transitions 
peaked at around 4.7 eV. As we have shown above these 
include the low-energy nominally dipole-forbidden, how- 
ever, phonon-assisted ii s (7r) — e g transition centered at 
around 1.7 eV, the weak dipole-allowed £211(71") — e g tran- 
sition centered at around 2.5 eV, and a more intensive 
dipole-allowed ti u (ir) — e g transition centered at around 
3.5 eV. All the p-d CT bands are characterized by a size- 
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able temperature dependence of the line-shape and spec- 
tral weight, both due to electron-lattice and spin-spin 
correlations, with a pronounced T-dependent SW trans- 
fer between dipole-allowed and dipole-forbidden bands. 
The situation becomes more involved, if to take into ac- 
count the interaction/mixing between p-d and d-d tran- 
sitions. 

Second, the analysis n ' 12 of the electron structure and 
energy spectrum of the d 3 -d 5 configuration to be a final 
one for the d(e 9 )-d(e ff ) charge transfer should be revisited 
in several points, in particular, as regards the assignment 
of the LS transitions. The fact is that the gap in be- 
tween the low-energy HS and LS d(e g )-d(e g ) CT transi- 
tions Ahs-ls = 10B + 5C«3 eV, i.e. markedly bigger 
than 2.3 eV suggested in Ref.HH. The authors have un- 
reasonably turned down the 4 E g ( 4 D) term of the d 5 con- 
figuration while the both 4 E g ( 4 G) and 4 E g ( 4 D) terms are 
almost equally involved in the d(e ff )-d(e ff ) charge trans- 
fer under consideration. As a result, both the assignment 
of the LS d-d CT transitions and quantitative estimates 
of such parameters as the Hund exchange integral J#, 
spectral weights SWls obtained in Ref.Ha seem to be 
questionable ones. 

Third, the authors address three bands at 2.0, 2.4, and 
2.7 eV to form a three-subband structure of the 2 eV d-d 
CT band. However, they cannot explain the origin of 
this splitting. Nevertheless, their quantitative analysis 
implies E\ = 2.0 eV (not 2.4 or 2.7 eV!) to be an energy 
of the HS d(e 9 )-d(e ff ) CT transition. 

From the other hand, the excellent experimental data 
by Kovaleva et a l 11 ^ 12 for LaMnC-3 are beyond all doubt. 
However, these support our assignment of p-d and d-d 
CT transitions in perovskite manganites and challenge 
a detailed analysis of the temperature effects with tak- 
ing into account a sizeable contribution of the electron- 
phonon coupling, comparable spin-dependent terms in 
the spectral weight for the both p-d and d-d CT tran- 
sitions, and coupling/mixing of the both types of CT 
transitions. 



B. Parent perovskite manganites: Mott-Hubbard 
or charge transfer insulators? 

Above we have provided an experimental evidence of a 
dual nature of charge gap states in perovskite manganites 
to be a superposition of inter-site d-d and one-center p-d 
CT transitions, a remarkable observation, which to the 
best of our knowledge has not been previously reported 
for the perovskite manganites. This implies a revisit of 
a standard ZSA classification scheme usually applied for 
these compounds. 

Zaanen, Sawatzky, and Allen (ZSA) 48 have pointed out 
that all the insulating 3d-electron systems can be classi- 
fied into two categories: the Mott-Hubbard and charge 
transfer insulators, depending on the nature of the gap. 
In the former the gap is formed by the transition be- 
tween the lower and upper Hubbard bands (U < A), 



while in the latter it is formed by the CT transitions 
between upper filled ligand and upper Hubbard bands 
(U > A). Actually in both cases we deal with CT transi- 
tions, however, in terms of a classification introduced by 
A. Moskvin et al. 22 i 23 these correspond to two-center (di- 
d2-) and one-center (pd-) CT excitons, respectively. Usu- 
ally one considers that the gap is of CT p-d type (U > A) 
for heavier transition metals (Cu, Ni, ...), while it is of 
d-d type (U < A) for the light transition metals (Ti, V, 
...). However, actual situation can be rather intricate. 
For instance, parent cuprates such as CuO, La2Cu04, 
Sr 2 Cu03, Sr 2 Cu02Cl2 are usually addressed to be Mott- 
Hubbard insulators, however, the detailed analysis of re- 
cent EELS data reveals its dual nature because of the 
gap described by a superposition of one- and two-center 
CT excitons 22 i 23 . Simple ZSA scheme needs in a revisit 
as it does not take into account a number of important 
issues. 

The parameter A p d appears in the simple generic pd- 
model describing the 3d oxides, with the hole Hamilto- 
nian 

H = J2 e ^d\ a d ia + £ P p]*Pjv + ^(tpddlaPj* + h.c.) 

i 3 ij 

+ Ud*^2 n ^n ld i + U P ^2 n 3P\ n ipi> ( 46 ) 
i 3 

where A p d = e p — e^, and all the terms have its usual 
form. This Hamiltonian describes fairly well the pd- 
hybridization effects. Depending on the sign and nu- 
merical value of the A p d parameter we come to a rather 
smooth redistribution of the hole density from the cations 
given large positive A p d to anions given large negative 
Apd. This effect, as it is stressed by Khomskii 4 ^ bears 
the quantitative character albeit the properties of sys- 
tems with positive and negative A p d can differ substan- 
tially. 

However, simple Hamiltonian (|4l)|) does not take into 
account a number of important effects, including the 
symmetry of different d-, p-states, crystalline field, polar- 
ization corrections, and electron-lattice coupling, partic- 
ularlly important in the case of a near degeneracy for the 
centers of gravity of the d- and p- manifolds. Electron 
correlations need in a more correct description, which 
must include a full set of the Slater, or Racah parameters. 
First of all, we should discriminate the (anti)bonding and 
nonbonding orbitals. For instance, in simple octahedra 
like Me06 we have two sets of 3d orbitals with e g (a) and 
t2 g (ft) symmetry, and seven sets of 2p orbitals with 
ai g ((r), h g (ir), ti u (a), t lu (ir), t 2u (^), e ff (cr), and t 2g (Tr) 
symmetry. Only two of them, e g and ti g hybridize with 
3<i orbitals having the same symmetry. From the other 
hand, namely the nonbonding O 2p orbitals seemingly 
form the first electron removal state in perovskites with 
octahedrally coordinated 3d ions—. Thus, simple clus- 
ter model points to a need in the realistic multi-band 
description of 3d oxides with a set of parameters like 
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A p< i— ■ The near degeneracy of the 3 c? orbitals with non- 
bonding oxygen states qualitatively differs from that for 
two hybridized orbitals of the same symmetry, and can 
result in a number of novel effects, in particular, gener- 
ated by enhanced electron-lattice coupling with local, or 
cooperative low-symmetry lattice deformations. In the 
case of the near-degeneracy of hybridizing 3d and 2p 
states the lattice deformation is governed by the active 
high-symmetry breathing modes. 

We should note that in many cases the charge trans- 
fer with the lowest energy corresponds to the dipolc- 
forbidden transition with a very weak spectral weight, 
that hampers its assignment. In general, we should make 
use of A p d symbol for the smallest of the O 2p- Me 3d 
charge transfer gaps, and > A pe j for the gap re- 

vealed in optical response due to strong dipole-allowed 
a — a and n — tt charge transfer transitions, respectively. 
Hence, one might say that the ZSA classification needs 
in detailization. 

As concerns the ZSA classification for perovskite man- 
ganites the situation seems to be far from being re- 
solved. The parent compound LaMnC>3 is sorted ei- 
ther into the charge transfer insulator—^ or the Mott- 
Hubbard insulator— & 11 ' 12 ! 17 , Interestingly, that in the 
study by Arima et alA LaMnC>3 is addressed to be lo- 
cated in the vicinity of the borderline where the optical 
gap changes from the Mott gap to the charge transfer 
gap that implies its dual nature. Our experimental data 
and theoretical analysis support this conjecture and ev- 
idence a dual nature of the dielectric gap in perovskite 
manganites RMn03, being formed by a superposition of 
p-d CT transitions and inter-site d-d CT transitions. In 
fact, the parent perovskite manganites RMn03 should 
rather be sorted neither into the CT insulator nor the 
Mott-Hubbard insulator in the ZSA scheme. 

At present there is no decisive experimental technique 
that provides the separation of one- and two-center CT 
excitons and distinction between Mott-Hubbard and CT 
insulators. It seems, one-dimensional 3d compounds are 
good candidates for such a study since the light with the 
electric field E perpendicular to the chain direction does 
excite only electron-hole pairs sitting on one MeO n clus- 
ter. On the other hand, for E parallel to the chain direc- 
tion, both types of excitons (one- and two-center's) can 
be observed. Indeed, the polarization-dependent angle- 
resolved EELS study of the ID cuprate Sr2Cu03 with the 
corner-shared CUO4 plaquettes has provided a unique op- 
portunity to separate both the one- and two-center CT 
excitons and reveal the two-peak nature of the CT gap 
with the presence of nearly degenerate two types of exci- 
tations—. In other words, similar to RMn03, insulating 
cuprates appear to be neither Mott-Hubbard, nor CT 
insulators, thus falling into intermediate region in which 
there are strong fluctuations between one- and two-center 
CT excitons. 

VI. CONCLUSIONS 

We have performed a compehensive theoretical and 
experimental study of optical response of parent per- 



ovskite manganites RMn03 . Starting with a simple clus- 
ter model approach we addressed both the one-center 
(p-d) and two-center (d-d) CT transitions, their polar- 
ization properties, the role played by structural parame- 
ters, orbital mixing, and spin degree of freedom. Optical 
complex dielectric function of single crystalline samples 
of RMn03 (R=La, Pr, Nd, Sm, Eu) has been measured 
by ellipsometric technique at room temperature in the 
spectral range from 1.0 to 5.0 eV for two light polar- 
izations: E || c and E 1 c. The comparative analysis 
of the spectral behavior of e\ and 62, particularly for 
E 1 c polarization, is believed to provide a more reliable 
assignment of spectral features. Our experimental data 
and theoretical analysis evidence a dual nature of the 
dielectric gap in these manganites, being formed by a su- 
perposition of forbidden ti g (Tr) —> e g (1.3 and 1.9 eV) or 
weak dipole allowed <2u(7r) — > e g (2.3 and 2.7 eV) p-d CT 
transitions and inter-site d(e ff )-de ff ) (2.0 eV in LaMn03 
and 2.2 eV in other manganites) HS CT transitions, a re- 
markable observation, which to the best of our knowledge 
has not been previously reported for these materials. In 
fact, the parent perovskite manganites RMn03 should 
rather be sorted neither into the CT insulator nor the 
Mott-Hubbard insulator in the Zaanen, Sawatzky, Allen 
scheme. We found an overall agreement between experi- 
mental spectra and theoretical predictions based on the 
theory of one-center p-d CT transitions and inter-site d- 
d CT transitions. Spectral features at 3.0 and 3.8 eV 
are assigned to a weak dipole-allowed fi u (7r) —¥ e g p-d 
CT transition of the it — a type with a clearly resolved 
low-symmetry splitting. Strong dipole-allowed p-d CT 
transitions ti u (ir) —} t2 g and t\ u (cr) — > e g of the 7r — 7r 
and a — a type, respectively, are observed as the most in- 
tensive bands with an unresolved structure at 4.7 and 5.3 
eV. Weak low-spin counterparts of the HS d(e ff )-de g ) CT 
transition are seemingly superimposed on these strong 
p-d CT bands. 

We argue that the both dipole-allowed d-d and p-d 
CT transtions are equally sensitive to the temperature- 
dependent spin correlations, only with a different relative 
change of the optical response through the onset of the 
magnetic order, relatively large for the former and rela- 
tively small for the latter. 

Finally one should note that some uncertainties and 
discrepancies observed in optical response of different 
i?Mn03 samples can be related with a charge transfer in- 
stability of perovskite manganites and a trend to a phase 
separation.— — 
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